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1. Introduction 

One of the most spectacular predictions of general relativity are black holes which should 
form when a large mass is concentrated in a sufficiently small volume. The idea of a 
mass-concentration which is so dense that even light would be trapped goes back to 
Laplace in the 18th century. Shortly after Einstein developed general relativity, Karl 
Schwarzschild discovered in 1916 a mathematical solution to the equations of the theory 
that describes such an object. It was only much later, with the work of physicists like 
Oppenheimer, Volkoff and Snyder in the 1930 's, that the scientific community began 
to think seriously about the possibility that such objects might actually exist in the 
Universe. It was shown that when a sufficiently massive star runs out of fuel, it is 
unable to support itself against its own gravitational attraction and it should collapse 
into a black hole. Starting with the 1960's and the 1970's, in the so-called Golden 
Era of black hole research, new interesting phenomena like the Hawking radiation and 
superradiance were discovered but for their rigorous mathematical description we have 
to wait until the 1990's and the beginning of the new century when the rigorous analysis 
of the propagation and of the scattering properties of classical and quantum fields on 
black hole space-times started to be developed. 

Whenever we attempt to analyze the scattering properties of fields in the more general 
framework of the Kerr-Newman black hole geometry, we arc faced with several difficulties 
which are not present in the picture of the Schwarzschild metric. First of all, the Kerr- 
Newman solution is only axially symmetric (cylindrical symmetry) since it possesses 
only two commuting Killing vector fields, namely the time coordinate vector field dt and 
the longitude coordinate vector field d,^. This implies that there is no decomposition 
in spin- weighted spherical harmonics. Moreover, another difficulty is due to fact that 
the Kerr-Newman space-time is not stationary. In particular it is impossible to find 
a Killing vector field which is time-like everywhere outside the black hole. In fact 
dt becomes space-like in the ergo-sphere, a toroidal region around the horizon. This 
implies that for field equations describing particles of integer spin (wave equation, Klein- 
Gordon, Maxwell) there exists no positive definite conserved energy. For field equations 
describing particles with half integer spin (Weyl, Dirac) we can find a conserved L2 norm 
with the usual interpretation of a conserved charge. Hence, the absence of stationarity 
in the Kerr-Newman metric is not really a difficulty for the scattering theory of classical 
Dirac fields. 

For the reasons mentioned above there are only few analytical studies of the propagation 
of fields outside a Kerr-Newman black hole. Time-dependent scattering for the Klein- 
Gordon equation in the Kerr framework has been developed by Hafner (2003), while 
Daude (2004) proved the existence and asymptotic completeness of wave operators, 
classical at the event horizon and DoUard-modified at infinity, for classical massive 
Dirac particles in the Kerr-Newman geometry. 

Our present work is the first of two papers devoted to develop a time-dependent 
scattering theory for massive Dirac particles outside the event horizon of a non-extreme 
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Kerr-Newman black hole. Before describing the contents of our work we should mention 
that the completeness of the Chandrasekhar ansatz has been proved for the first time 
in Finster et al. (2003). There, they derived an integral representation for the Dirac 
propagator in terms of the solutions of the radial and angular ODEs arising by means 
of the Chandrasekhar separation of variables which in turn they used to read out the 
completeness of the ansatz. In the present work we proceed in the opposite way. The 
main advantage is that the derivation of the propagator will follow quite immediately 
after we have proved the completeness of the Chandrasekhar ansatz. As a by-product, 
we also obtain an integral representation for the Dirac propagator in the Minkowski 
space-time in oblate spheroidal coordinates. Our integral representations together with 
the estimates for the asymptotic behaviour of the solutions of the radial problem will 
allow us in the second paper to give explicit analytical expressions for wave operators 
classical at the event horizon and DoUard-modified at infinity. Asymptotically away 
from the black hole it will turn out to be of advantage to consider first classical 
wave operators since their computation will lead us to an analytical expression for 
the time-dependent logarithmic phase shift needed to construct the DoUard-modified 
wave operators. After implementation of this phase shift in the free dynamics we will 
evaluate the Dollard-modified wave operators and obtain an integral representation for 
them. To our knowledge this will be the first analytical result since all previous works 
in this direction (see Hafner (2003) and Daude (2004)) which are based on the Mourre 
theory (Mourre (1981)), treat only the problem of the existence of the above mentioned 
wave operators. 

The present article is organised as follows: In Section 1 we motivate our interest in the 
derivation of an integral representation for the Dirac propagator in the Kerr-Newman 
metric. In Section ITT] we give a short review of the Dirac equation in the Kerr-Newman 
geometry. In Section 1?!^ we introduce the so-called Chandrasekhar ansatz. In Sectional 
we compute the scalar product with respect to which the Dirac Hamiltonian obtained 
from ()2.5p is formally self-adjoint. In Section 0] we bring the modified but equivalent 
Dirac equation ()2.15j) into the form of a Schrodinger equation with Hamiltonian h and 
we show that h defined on the set of the smooth functions with compact support 
is essentially self-adjoint. In Section El the main results are Theorem 15.21 where we 
prove the completeness of the Chandrasekhar ansatz and Theorem 15.41 which gives the 
integral representation for the Dirac propagator in a Kerr-Newman black hole manifold. 
Lemma lFTTl and l(j.2l in SectionlHldescribe the asymptotic behaviour of the radial functions 
satisfying the system of first order ODEs ()2.36|) introduced in Section 12.21 Finally, in 
the last section we also derive an integral representation for the Dirac equation in the 
Minkowski space-time in oblate spheroidal coordinates. 
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2. Preliminaries 



2.1. The Dirac equation in the Newman-Penrose formalism 

According to Penrose and Rindler in terms of two-component spinors (0"^, x"^ ) the Dirac 
equation coupled to a general gravitational field and a 4- vector field V is given in Planck 
units h = c = G = lhY 



nie 



(V 



teVi 



ITLf. 



[2.1) 



where V^^' is the symbol for covariant differentiation, e is the charge or coupling 
constant of the Dirac particle to the vector field V and mg is the particle mass. Notice 
that the factor 2" 3 appearing in the above equations is due to the fact that the Pauli 
matrices as defined in the Newman-Penrose formalism differ from their usual definitions 
by the factor \/2. The next step is to bring ()2.H1 into a system of partial differential 
equations with respect to a coordinate basis. In the framework of Newman and Penrose 
this is achieved by choosing a null tetrad, i.e. a set of four vector fields (1, n, m, m) 
where (1, n) is a pair of real vectors and (m, m) a pair of complex conjugate vectors. 
Moreover, they satisfy the orthogonality conditions l-m = l- m = n- m = n- m = 0, 
the requirements that the vectors be null, i.e. l-l = n- n = m- m = m- m = and 
finally, the normalization conditions 1 ■ n = 1 and m • m = —1. The covariant derivative 
is decomposed into directional covariant derivatives along the frame vectors which we 
denote by 

D = rVa, A = n'^Va, 5 = m'^Va, 5 = rTVa. (2.2) 

The spin coefficients can be written in terms of the Ricci rotation-coefficients (e.g. 
Chandrasekhar (1992)) by defining the frame vectors and their dual 1-forms as follows 
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where enclosure in parenthesis distinguishes the tetrad indices from the tensor indices. 
The Ricci rotation-coefficients are expressed by 
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\a)(b)(c) = e(b)i;j 



where g is the metric tensor. The spin coefficients can be written by means of the 
7(a) (6) (c) as follows 
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7(3)(1)(1), 
7(3)(1)(3), 
7(2)(4)(4), 
7(2)(4)(2), 
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T(3)(l)(4), 
■ 7(2)(4)(3), 

7(3)(1)(2), 
^ 7(2)(4)(1), 



[7(2)(1)(1) +7(3)(4)(1)] 



7=2 h2)(l)(2) + 7(3)(4)(2)] 
a = - [7(2)(1)(4) + 7(3)(4)(4)] 
/3 = - [7{2)(1)(3) + 7(3)(4)(3)] 



(2.3a) 
(2.36) 
(2.3c) 
(2.3d) 



Let us recall that to any tetrad (1, n, m, m) we can associate a unitary spin-frame 
{o^,t^), defined uniquely up to an overall sign factor by the requirements that (for 
details see Penrose and Rindler) 

Denoting by 0o and (pi the components of in {o^, t^) and by Xo' and Xi' the 
components of x^ in {o^ > 'i^^ ) 

00 = 0ao^, 01 = (pAt"^, Xo' = Xa'O^ , Xi' = Xa''^^ ' 
the Dirac equation ()2.1|) takes the form 



{D + e- 


P- 


- ierVa)Xi' - 


- (5 + TT — a - 


- ienfVa)xQ' 




(2.4a) 
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- iefn''Va)Xi' 




(2.46) 


(A + /i- 


7 


- ien"K)0o - 


-i6 + f3-T- 


- zem"K)0i = 


me 


(2.4c) 
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Let us introduce the following symbols 

a+ = -iV2{D + e - p - ietVa), (3+ = i^/2{6 + ¥ - a - iem^Va), 
/?_ = iV2(5 + P-T- ieWVa), a_ = -iv^(A + /I - 7 - ien^K), 
5_ = -iV2(A + p - 7 - ien^K), ^+ = «V2('5 + /? - r - iem^'K), 
^_ = ^^2(5 + TT - a - iem"\4), 5+ = -iV2{D + e- p- ietVa) 
and let us define the spinor \I' = (Fi, F2,Gi, 6*2)^ with components 

i^i = -00, = 01, 6*1 = iXi', G'2 = iXo'- 
Then the system of equations (2.4a)-(2.4c?) is simply 

/ —rrie q;+ /3+ \ 
—me /?_ a_ 
5- — /3+ —me 

y — 5+ —me y 

In what follows, we consider the Dirac equation in the Kerr-Newman geometry, i.e. in 
the presence of a rotating charged black hole. In Boyer-Lindquist coordinates {t, r, -d, ip) 
with r>0, 0<'i9<7r, 0<(/)<27r the Kerr-Newman metric is given by (e.g. Wald) 

, , A - sin^ ?9 , , 2a sin^ ^r"^ + a^ - A) , , 
ds^ = dt^ H ^ -dtdip 

_ ^dr^ - T.d§^ - (r^ + a'f sin^ ^(2.5) 

with 

E := S(r, e)=r'^ + c? cos^ A := A(r) = - 2Mr + a^ + 



^ = 0. (2.5) 
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and 

E := E(r, -d) = 1 — a^7^(r) sin^ 'd, 7(r) : = 



+ 



where M, a and Q are the mass, the angular momentum per unit mass and the charge 
of the black hole, respectively. Here we will always consider a 7^ 0. In the non-extreme 
case > + function A has two distinct zeros, namely, 

the first one corresponding to the Cauchy horizon and the second to the event horizon. 
Moreover, A > for r > ri. In the extreme case = + ^]-^g Cauchy horizon and 
the event horizon coincide since A has a double root at = M. 

Lemma 2.1 Let > a"^ + Q"^. For every r > ri > and i} G [0, vr] it results S > 0. 
Proof. Let us write S as follows 

E=-4^^, p{r,{}) = {r' + ay -a' Asm' 

A short computation gives 

p(r, '&) = + a'{l + cos^ ^9)r^ + 2MaV sin^ + a^(a^ cos^ - Q'^ sin^ -d) 
where we used the definition of A. Since r > ri > M, we have 

p{r, d)> + a'M'il + cos^ + 2M''a' sin^ § + 0^(02 cos^ ^ - Q'^ sin^ ^9), 
= + 2a' M'^ + cos^ ?9 + a'{M' - Q^) sin=^ §. 

Recall that in the non-extreme case M' > Q', while in the extreme case M' — Q' = a'. 
Hence, the last line in the above expression is positive for every G [0, tt]. □ 

Finally, notice that by setting M = Q = in ()2.6|) the Kerr-Newman metric goes over 
to the Minkowski metric in oblate spheroidal coordinates (OSC), namely 

ds' = dt' - ^dr' - Udd' - A sin^ dd^' , A = r' + a' . (2.7) 
A 

In fact, by means of the coordinate transformation 

X = ^/r"^ + a' sin ^ cos y = \/ r"^ -\- a' sin "d sin z = rcos^ (2-8) 

fl2.7|l can be reduced to the Minkowski metric in cartesian coordinates. In the case 
M = Q = the function A has no real zeros. Thus there are no coordinate singularities. 
However the metric has a curvature singularity at S = 0, i.e. for r = and d = ^. 
The surfaces of constant r are confocal ellipsoids which degenerate at r = to the disc 
z = and x' + y' < a'. The surfaces of constant d are represented by hyperboloids 
of one sheet, confocal to the above mentioned ellipsoids and ^ = f corresponds to the 
boundary of the disc at x' + y' = a'. As a consequence the curvature singularity occurs 
on the boundary of the disc, i.e. on the circle x' + y' = a' and z = 0. Moreover, 
there is no reason to restrict r to be positive since the space-time can be analytically 
continued through the disc to another flat region with r < 0. For further details we 
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refer to Flammer (1957) and to the original work of Boyer and Lindquist (1967). 
The Dirac equation in the Kerr-Newman geometry was computed and separated by 
Page (1976) with the help of the Kinnersley tetrad (see Kinnersley (1969)). We recall 
that such a tetrad is naturally inherited from the Petrov type D structure of the Kerr- 
Newman metric (e.g. O'Neill (1995)). In view of the separation of the Dirac equation 
we choose to work with the Carter tetrad (Carter (1987)) 

rV„ =^(r^-~-Al + l) (2.9a) 

V2AS V dt dr dip J 

m'^Va = , ^ { iq^^ - {a' - q')^ - ^-^] (2.9c) 

/2(a2-g2)s V dt ^ 'dq 3^ J ^ ' 



where the dependency of the new variables t, q and (p on the Boyer-Lindquist variables 
is given by the relations 

t = t — a(p, g = acos'i9, ip = —. 

a 

In the above defined variables the electromagnetic potential is simply given by 

Va = ^{lMq'), ^ = r' + q' 

which is stationary and axially symmetric. Using the above tetrad and ()2.3aj) -f2.3(j) the 
spin coefficients are computed to be 

K = (j = A = z/ = 0, p = /i, 7 = e, 71" = — T, [5 = —a, 

i{a'^ — q^)r — (r^ + a^)q 

a = 

2v/2(a2 -g2)s 

with p = —{r — iq)~^. Notice that the vanishing of k and a implies that the congruence 
formed by the vector 1 is geodesic and shear-free, respectively. As a consequence the 
Riemann tensor is of type 11 according to the Goldberg-Sachs theorem (see Goldberg 
and Sachs (1962)). 

We are now ready to give a more explicit form to the Dirac equation ()2.5|1 . Indeed, by 
introducing the following operators 



V - — - 

^ dr A 



(2.10) 



= ^ + 1 cot^ T . (asin^l + csc^^^ , (2.11) 
the entries of the matrix in ()2.5|) are computed to be 

a± = ±i^{V± + f{r,{})), /5± = ^(£i + <^(r,^)), (2.12a) 



5± = ±^A/^(P± + /(r,^9)), []^ = -^{C^ + g{r,d)) (2.126) 
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1 /r — M 1 \ , ia sin {} 

/(^' ^) = o — ^ + zrr—zzzTo ' 9i^^ ^) = -TT?:-— — -TV- (2-13) 



with 

2\ A r + iacos'i?/' ' 2(r + iacos'i?) 

r/ie Chandrasekhar separation ansatz 

Carter and McLenaghan (1979) found that in the Kerr geometry and in the fiat-space 
limit there exists a generalized total angular momentum operator of the form 

^ = ^757^^ (/;^. - ^7 , D, = V,-teA,, (2.14) 

commuting with the Dirac operator under the condition that is an antisymmetric 
tensor satisfying the Penrose-Floyd equation f^{u;p) = (see Penrose and Floyd (1973)). 
In particular Carter and McLenaghan also showed that ()2.14|1 can be interpreted in the 
Minkowski space as the square root of the ordinary total squared angular momentum 
Casimir operator of the rotation group. In this section we recall how to construct from 
the Chandrasekhar ansatz (see Chandrasekhar (1976)) for the Dirac equation in the 
Kerr-Newman metric the symmetry operator J for the formal Dirac operator Oo defined 
in ()2.5j) and we compute it explicitly. As pointed out in Carter and McLenaghan the 
existence of such an operator underlaying the Chandrasekhar separability of the Dirac 
equation is due to the presence of an appropriate Killing spinor field on the space-time 
under consideration. For the concept of Killing spinor we refer to the original work of 
Penrose and Walker (1970). 

We begin by replacing the Dirac equation ()2.5|) by a modified but equivalent equation 

= {TS-^OdS)^^ = (2.15) 

where tjj = S^^"^ = {Fi, F2, Gi, 6*2)"^ and F and S are non singular 4x4 matrices, whose 
elements may depend on the variables r and "(9. The next lemma shows how to find non 
singular matrices 5* and F such that the formal operator W decomposes into the sum of 
two formal operators, the first one containing only derivatives respect to the variables 
t, r and (f and the second one involving only derivatives respect to t, and ip. 

Lemma 2.2 Let r he positive with r > ri. Moreover, let S and F be non singular 
matrices defined as follows 

S = A-3 diag ( , ^ ^ , ^ ^ ) (2.16) 

V V — cos yr — iacos^ yr + ia cos-d yr + iacos-dj 

with det(5') = (SA)^-'^ and 

V = —i diag(r + ia cos -(9, — (r + ia cos -(9) , — (r — ia cos 'd),r — ia cos -d) 
with det(F) = S^. Then 

W = >V(t,,,^) + >V(t,^,^) (2.17) 



where 



( 



(t,r,(/p) 
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—airte cos 

















arrie cos 
















arrie cos 



















arrie cos 



(2.18) 



(2.19) 



with V± and C± defined by \2.1U\) and \2.11]) in Section WH 

Proof. Let us introduce an invertible matrix 

S := diag(/;,(r,79), A(r,'i9),cr(r, ^9), 7(r, -i?)) 

with functions h, A, o", 7 at least C"'^((ri, +00) x [0, vr]). Then the equation OdStP = 
gives rise to the system 

- rriehFi + a+{(jGi) + P+{-fG2) = 0, (2.20a) 

- meAFs + /3_(aGi) + a_(7G2) = 0, (2.206) 

- rrieaGi + a_{hFi) - p+{AF2) = 0, (2.20c) 

- mejG2 - ]3-{hFi) + 5+(AF2) = 0. (2.20(i) 

If we substitute ()2.12a|l -(2.126). we obtain 



a+(o"Gi) = zi 
Let a and 7 be 



(97 



(2.21a) 
(2.216) 



5(^) 



Aiy/r + ia cos?? 



, 7(r,79) 



7(r) 



V^H-iocosl? 

with the constraints 5('i9) 7^ for every G [0, tt] and 7(r) 7^ for every r G (ri, +00). 
Then (j2.21a] -M simphfy to 

^ /a" ^ i ^ 

a+{(rGi) = id —aV+Gi, P+{'yG2) = ^=7£+G'2- 



The substitution of (I2.12al) -r2.126) in (lOIl gives 



"-(7^2) 



/?-(^Gi) 



72^- + 1^ + /7 1 ^2, 



(2.22a) 
(2.226) 
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It is easy to verify that (|2.22aj) -f2.226) reduce to 



by choosing 7(r) = C4A 4 with C4 G C\{0} and a{^) = C3 E C\{0}. Analogously we 
find that 



-, A(r,^) 



Aiy^r — ia cos'd ' A4 + ia cos"!? 

for some Ci, C2 G C\{0}. Let us define a non singular 4x4 matrix 

T = —i diag(r + ia cos {},—{r + ia cos {}), —{r — ia cos 'd),r — ia cos ^9) . 

After a tedious computation we find that 



w = rs-^os 



( 



\ 



imer—anie cos i 


C4 







-imer+arrie cos i 



/AX'4 



C2 



C2 





-irrisr—ame. cos i9 

imer+arrie cos 1? J 



If we require that ci = C2 = C3 = C4 =: c, we finally have 
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im^r—ame cos i 


— 





-imer+ame cos i 



— £_ 



-zmer— ame cos i 






imer+arde cos •& 



Without loss of generality we can set c = 1 and this completes the proof. □ 

Remark. It should be noted that the regular and time-independent transformation S 
(eq. (2.1) in Finster et al . (2003)) is in contrast with our ()2.16|) . 

Notice that Lemma [2.21 is still valid also in the extreme case where A = (r — M)^. Since 
by setting M = Q = (see discussion in Section ITTj) the Kerr- Newman metric becomes 
the Minkowski metric in oblate spheroidal coordinates, we obtain from Lemma 12.21 an 
analogous result for the Dirac equation in OSC which we state in the following corollary. 



Corollary 2.3 Let TZ denote the circle TZ = {{x,y,z) G 
Moreover, in M.^\R let us define the non singular matrices 

S = A 4, diag 



a?, z 



1 



0}. 



^\Jr — ia cos ' \Jr — ia cos??' Vr + ia cos"!?' Vr + ia cos?? 
with A = + a^ and F as in Lemma \2. M Then 



y^,OSC, ^ ^iosc, + ^^^^^^^^ 
/ im^r 



^{OSC) 



—irrier 
AP_ 

\/ap. 



'AV+ 


—irrier 
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AV_ 


irrier ) 



(2.23) 



(2.24) 
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where Wu -d is given by \2.19^) and 



^ (9 f d ad 

It is interesting to observe that the angular operator W(t,^,^) remains unchanged for the 
Dirac equation in the Kerr-Newman metric and in the Minkowski space-time in oblate 
spheroidal coordinates. We compute now the commutation relations which we will need 
later to verify that in the Kerr-Newman metric the operator J given by (|2.14j) is a 
symmetry operator for Od defined in ()2.5|) . It is straightforward to check that 

[W(t,.,^), W(t,^,^)] = 0, [W(t,.,^), W] = >V] = 0. 

Furthermore, the matrix F introduced in Lemma f2.2l splits into the sum 

r = F(.r) + F(^) 

with F(r) and F(^) satisfying the commutator relations 

[r(r),F(^)] =0, \^{r),y^{t,'d,Lp)\ = [^(^),yV(t,r,ip)] =0. 

In what follows we restrict our attention to stationary solutions of the Dirac equation. 
In particular since the Kerr-Newman metric is axially symmetric, it is natural to make 
the following ansatz for the spinors ip entering in ()2.15j] 

^{t, r, ^, if) = e*"*e'('^+^)^^(r, t?) (2.28) 

where uu and G Z denote the energy and the azimuthal quantum number of the particle, 
respectively. Let us substitute ()2.28|1 into ()2.15|1 . Then ip{r,^) satisfies the equation 

+ W(^)) ^ = (2.29) 

where 



(2.25) 
(2.26) 
(2.27) 



(r) 



inier VAX'+ 

-imer VaS- 

VA23_ -inier 

\^'D+ imer 

^ —am.ecos'd 

anie cos-d — 

C+ —a»ne cos 1? 

—iZ— arrie cos 



(2.30a) 



(2.306) 



with 



K{r) = uj{r^ + a^)-eQr+ ik + -]a, (2.31a 



^ _ 9_ ^. K{r) 
^ dr A ' 

C^= + cot^±Q{^), Q{^) =aLJsm^+ (k + -]csc^. (2.316) 
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If we set ifjir^'d) = {fi{r,'d), f2{r,'d),gi{r,'d),g2{r,{}))'^, equation (j2.29j) gives rise to the 
following system of first order linear partial differential equations 

VAP+5fi + inierfi + {£,+§2 - am^ cos^9/i) = 0, 
\fKV_g2 - im^rj2 - (l^-gi - arrie cost9/2) = 0, 
\^V_fi — imeVgi + (£+/2 — arrie cosi!}gi) = 0, 
\^V+f2 + imerg2 — {C-fi — anie co?>dg2) = 0. 
Let us now define 

A(r,t9) = 7iW5i(^), /2(r,t9) = a2(r)r2(^), 

gi{r,§) = ai(r)/3i(?9), g2{r,^) = e2(r)/i2(^9). 
It can be easily seen that the separability conditions for the above system are 

AW = 6i{^), e2(r) = 7iW, /X2(^9) = ^2(^9), ai(r) = a2{r). 
Following Chandrasekhar (1976), we set 

7i(r) = i?_(r), SM=S-m, a2(r) = i?+(r), r2(^) = 



and we get for tp{r,i)) the so called Chandrasekhar ansatz 

/i(r,t9) = i?_(r)5_(7?), /2(r,^) = R+{r)S+{^), (2.32a) 

g,ir,^) = R+{r)S4^), g2{r,^) = R-{r)S^{^), (2.326) 

by means of which the above system of partial differential equations becomes 

{y/AV+R+ + imerR^)S^ + (£+5+ - am^ cos^S^)R^ = 0, (2.33a) 

(VAP_i?_ - im^rR+)S+ - (£_5_ - am^ cos 'dS+)R+ = 0, (2.33&) 

{^/KV_R_ - im^rR+)S^ + (£+S+ - ame cos ^9S_)i?+ = 0, (2.33c) 

{y/AV+R+ + imerR-)S+ - (£_5_ - ame cos t95+)i?„ = 0. (2.33d) 

Introducing four separation constants Ai, . . . , A4 as follows 

£+5+ — arrie cos ■i}S^ = — XiS^, (2.34a) 

C^S^ — ameCos'&S+ = + A25'+, (2.346) 

C-^-S+ — ameCos')!}S- = — X^S-, (2.34c) 

C^S^ — ameCos'&S+ = + A4S'+, (2.34c?) 

we obtain from ()2.33a|) -f2.336?) 

\/~Kv+R+ + irUerR^ = XiR_, (2.35a) 

\^V+R+ + irrierR. = X^R., (2.356) 

VAP_i?_ - imerR+ = X2R+, (2.35c) 

VAP_i?_ - imerR+ = XsR+. (2.35d) 



Clearly, the systems of equations ()2.34a|l -(2.34d) and ()2.35a|l -(2.35d) will be consistent 
if Ai = A2 = A3 = A4 =: A. Notice that the decoupled equations give rise to the 
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following two systems of linear first order differential equations for the radial and angular 
components of the spinor ip 



0. (2.37) 




irrier — A 

—C- X + arUe cos 6 
X — anip cos 6 C+ 




For a = the components S± of the angular eigenfunctions can be expressed in terms 
of spin-weighted spherical harmonics (see for details Newman and Penrose (1966) and 
Goldberg et al . (1967)) whereas in the more general case a 7^ it can be shown that 
they satisfy a generalized Heun equation (see Batic et al . (2005)). 
Starting from the systems of equations (|2.34ap -f2.34c?) and (|2.35aj) -f2.35d). we can 
construct the symmetry operator J for the Dirac operator On defined in fl2.5|) . First 
of all we set Ai = A2 = A3 = A4 =: A in (jZMi|)-(2.34(i) and (IZSSi|)-(2.35c?). 
After multiplication of ()2.34a|) and ()2.35aj) by — ri?_ and iacosdS^, respectively, and 
summation of the resulting equations we obtain 

— iap cos ^yfKv j^gi + rpC+g2 = A/i (2.38) 

with p = — (r + ia cos Proceeding analogously we obtain from the remaining 
equations 

— iap cos d\/~K'D_g2 — rpC^gi = A/2, (2.39a) 
+ zap cos ?9VAP_/i + rp£+/2 = Xgi, (2.396) 
+ iap cos i9\^V^f 2 — rpC-fi = Xg2, (2.39c) 

with p = — (r — iacosi!})^^. By means of ()2.38p and ()2.39al -c) we can obtain an explicit 
expression for the matrix operator J, namely 

f -ia^ cos ^VAV+ rpC+ \ 

—rpC- — iap cos 'd\/~KV_ 

iap cos ■&^/AV_ t"P^+ 

y —rpC_ iap cos'&^/AV^ J 



J 



According to Carter and McLenaghan J can be written in a more compact form as 
follows 

j = r-i(r(^)W(t,.,^) - r^r)WitM). (2.40) 

Starting from J the next lemma shows how to construct a new operator J that will 
commute with the Dirac operator ()2.5|) . 

Lemma 2.4 Let J be defined as in \2.40i )- Then in the Kerr- Newman metric the 
operator J = S JS^^ with S defined as in Lemma [2.^ is a symmetry matrix operator for 
the formal Dirac operator Oo introduced in \2.3) . i.e. 

[Od, J] =0. 
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Proof. From (|2.15|) we have Od = ST ^WS ^ acting on the spinor Let us consider 
the operator J = SJS^^ with J given by ()2.40j) . Then we find that 

OdJ = sr-iwr-i(r(^)W(t,,,^) - r(,)W(t,,,^))5-i 
jOd = 5r-i(r(^)>V(t,,,^) -r(,)W(t,^,^))r-^ws-i 

and therefore [Od, J] = ST-^PS-^ with 

Making use of the decomposition W = W(r) + W(^) and of the commutation relations 
(IT^ . we obtain 

P = >V(i,^,^)(r(,)r-i + r-ir(^))W(t,,,^) - W(t,,,^)(r-ir(,) + r(^)r-i)W(t,^,^). 

Finally, with the help of ()2.25|) it results that 

P = W(t,tf,^)W(t,,,^) - W(t,,,^)W(t,^,^) = 
and this completes the proof. □ 

3. The physical scalar product 

The Dirac equation ()2.5|) possesses a conserved current (e.g. Wald) on general curved 
space-times. As a consequence the total charge outside the black hole 

Co = C{to) = [ {cI^aJa' + XaXa')T^^ dV 

is constant throughout time where denotes the future oriented normal vector field 
to a Cauchy surface Sj^, normalized so that T^T" = 2 and dV is the volume form 
induced on by the metric g. Here, a denotes an abstract tensor index which has to 
be understood as an unprimed spinor index A and a primed spinor index A' clumped 
together, i.e. a = AA' . We recall that for a given null tetrad the vector field + 
is time-like future oriented since it is the sum of two future-oriented null vectors. This 
allows us to associate to any null tetrad a preferred time-like future-directed vector field 
1°- + n"'. Thus, it is natural to impose = + n°- which becomes in the Carter tetrad 
and in Boyer-Lindquist coordinates 

2ir' + a') fd , a d 



a/2AS \dt r'^ + a?dip^ 
Taking into account that for a generic null tetrad it was shown by Hafner and Nicolas 
(2004) that 



T 



AA 



and employing ()2.9aj) -f2.9c) a simple computation gives that the matrix T^^ is the 
identity. Therefore the total charge outside the black hole is 

Co= I {\W + \<Pi? + \Xo'? + \Xi'?)dV 
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which is manifestly a positive quantity. By means of Theorem 13.3 in Oloff the volume 
form dV induced on J^to by the Kerr-Newman metric is computed to be 

/»+oo /»1 /'27r 

Co= / dr d{cos'd) / dip{r^ + a^) 

Jri J-1 Jo 

where the spinor \1/ has been defined in Section ITTl Let us bring the Dirac equation (j2.5|) 
into the form of a Schrodinger equation idt"^ = Hd^! . Then, the formal Hamiltonian 
operator Hd acting on the spinor \1/ on the hypersurface is formally self-adjoint with 
respect to the positive scalar product 

dr j d{cos^) j c/v;(r2 + a2)y — (3.1) 

If we introduce the tortoise coordinate m G M defined by du/dr = (r^ + a^)/A, the scalar 
product p.ip can be written in the more compact form 

/+00 nl i>2n ^_ 

du d{cos^) rfy^VSSA^^. (3.2) 
-oo J-1 Jo 

Finally, by trasforming the spinors according to \I' = (SSA)3vl/ and taking into account 
the corresponding changes in the Hamiltonian Ho due to this transformation, we could 
even work in a standard (i.e. not weighted) L2 space. 

4. Reduction to the Schrodinger form and self-adjointness of the Dirac 
Hamiltonian 

In order to bring the matrix equation Wtp = with W given by ()2.17j) into the form of 
a Schrodinger equation idt^ = hijj we bring the time derivatives in ()2.17|) to the l.h.s. 
of the equation and find that 

- zldti^ = (W(,,^) + W(^,^))V^ (4.1) 



with 



/ -asm^ \ 



1 



— asini? —i 



/A 



-asint? 
y -asin^9 ^ / 



and W(r.,(^), VV(^,^) formally given by ()2.18p and ()2.19|) . respectively, with V± and £-1- 
now replaced by 

^± = ^ T ^ ^a- - ^eQr^ , = - + - cot ^ ^ . esc ^- . 
Lemma 4.1 det(T) 7^ for every r > ri > and ^ E [0, tt]. 

Proof. A short computation gives det(T) = p^(r, '(9)/A^ with p{r,d) = (r^ + a^)^S. 
Since Lemma (2.11 implies that in the non-extreme case (as well as in the extreme case) 
E can never be zero, it follows that T is invertible. □ 
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Let u be the tortoise coordinate introduced in Section |21 According to Lemma 14.11 we 
can multiply 1)4.11) on both sides by — to obtain 

h = ho + V{u,^) (4.2) 

where the formal differential expression Hq is given by 

\ / 



ho = A{u,'d) 

















£+ 




























-s 



+ 



and 



1 
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a'-y{u) sin-i? 






-M+ 




























M 














/J 



(4.3) 




















where ct2 is a Pauli matrix, p 


= — (r — 


ia cos 


?9)-i 



eQr 



(4.4) 
(4.5) 



^ . d ia d 

^ du r'^ + dip ' 



M± = i'y{u)L 



satisfying £± = —£± and M.j^ = —M.^. In the definition of S± we made use of the 
tortoise coordinate u e (— cxd,oo), that we introduced in Sectional Furthermore, the 
matrix contained in the square brackets in ()4.4|) is hermitian. 

Remark. It should be noted that the product of the matrices B and C entering in the 
equation (2.11) in Finster et al . (2003), more precisely 



BC 



is in contrast to our result ()4.4|) . 

In view of fl3.2|) in Section El it follows that the formal Hamilton operator h is formally 
self-adjoint with respect to the positive scalar product 

"+00 




du 



d^cos'd) / dip 



2tt 



(f) (4.6) 
Sijj with the transformation 



where the spinors has been transformed according to \I' 
matrix 5* given in Lemma [2.21 and E defined in Section ITTl 

Remark. As in Finster et al . (2003) we made use of the Carter tetrad and after a 
regular and time-independent transformation given in Lemma 12.21 we found that the 
Dirac equation ()2.17j) coincide with their equation (2.1). However the scalar product 
expressed by (2.15) and (2.16) in Finster et al . (2003) does not agree with ours ()4.6p . 
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For simplicity in notation we will write instead of when there is no risk of 
confusion. In what follows we consider the Hilbert space Ti. consisting of wave functions 
ip : Q := X S"^ — > together with the positive scalar product ()4.fjj) . We study now 
the formal differential expression Hq as an operator on 

To this purpose we introduce the domains D{Hq) consisting of all G L2(fi)^ such that 
Hqi/j G L2{^y in the sense of distributions and -D(T) := C^{Qy where HqiP = hoip for 
all e D{Ho) and f^j = hoi^ for all e D{f). 

Lemma 4.2 T is essentially self-adjoint and T = Hq. 

Proof. We show first that T is symmetric. Clearly for all ip, (p & D{T) the operator T 
is formally self-adjoint with respect to the scalar product ()4.6|) . i.e. 

{tm = (V^|T0). (4.7) 

Since T is densely defined and fl4.7|l holds for every ip, (p E D{T), it follows that T 
satisfies the definition of a symmetric operator as given in Birman and Solomjak. The 
proof of the essentially self-adjointness is the same as that for part (c) of Theorem IX.27 
in Reed and Simon vol II. □ 

Notice that, since T is essentially self-adjoint, then there is uniquely associated to T a 
self-adjoint operator T = T** = T* = Hq and we call it the unperturbed Hamiltonian. 
Let S be the minimal operator constructed from Hq + V, i.e. S = T-\-V = T + Vo + Vi 
where V, Vo and Vi denote the maximal multiplication operators by V, Vq and Vi with 

V{u,^) = Vo{u,^) + Vi{u,^), 

Theorem 4.3 S is essentially self-adjoint and S = H . 

Proof. A simple calculation shows that Vq, Vi G Loo(f^) since 

\eQ\ 

llFolloo := sup{|yo(a;)|} < 4^^, \\Vi\\oo := sup { 1 14 (x) | } < 4me. 

Hence, Vq and Vi are bounded and S and T have the same domain C^{VlY. Let us 
consider S as obtained from T by adding the perturbing term V. In order to show 
that S is e.s.a., it suffices to prove t_hat V is relatively bounded with respect to T with 
relative bound < 1. Let H denotes 5, we show that 

H = Ho + V, D{H) = D{Ho) C D(y). 

In view of (5.12) in Kato (Ch.V §5.3) since Viip G L2(fi)^, we have 

llVi^lU, < III^iIIooII^IIl, <c||Fi|U(a-i/l/7o^|U,+a^/l^|UJ 
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with some constant c > and arbitrary a > 0. On the other hand also Voip ^ -^2(^)^ 
and llVoV'lUa < IIK)||oo||^||l2- Hence, we have D(V) D D{Ho) D D{f) and 

IIV^IU, <a||^|U,+6||i/o^|U, 

with b = ca""*^/^!! Villoo and a = ca'^/^H Vi||oo + ||Vo||oo- Since we can choose a arbitrary 
large, the above inequality shows that V is i^o-bounded with relative bound 0. Hence, 
Theorem 4.1 in Kato (Ch.V §4.1) implies that S = T + V is essentially self-adjoint. 
Finally, since Hq + V D 5, it follows that Hq + V = S. Thus, the perturbed operator H 
has the same domain as the unperturbed operator Hq. □ 

Remark. It should be noted that in Finster et al . (2003) the essentially self-adjointness 
of H is stated without proof. 



5. The Dirac propagator in the Kerr-Newman metric 

We begin with some preliminary observations. Energy, generalized squared total 
angular momentum and the z-component of the total angular momentum form a set 
of commuting observables {if, J^, Jz}- Moreover, the angular system ()2.37|1 can be 
brought into the so-called Dirac form (see Batic et al . (2005)), namely 

, , ^ / l\ dS ( —aiTLe cos -i? — -^T^^ — acusmd \ ^ , ^ , ^ , 
US := I T n 7i + k+i KS = (5.1) 

y ~1 J dv y ___| _ a^jsin?? anie cos i!) J 

with S{i}) = {S^ (??), S'+('i9))^ and i} G (0, it). For the theory of Dirac systems of ordinary 
differential equations we refer to Weidmann (1987). In L2(0,7r)^ the angular operator 
U with domain D(U) = C^(0,7r)^ is essentially self-adjoint, its spectrum is discrete and 
it consists of simple eigenvalues, i.e. Xj < Xj+i for every j G Z\{0}. Moreover, the 
eigenvalues depend smoothly on u. For details see Finster et al . (2000) and Batic et 
al . (2005). Furthermore, the functions e^^^^^^'^ are eigenfunctions of the z-component 
of the total angular momentum operator Jz with eigenvalues ~ + |) with A; G Z. 
Hence, in the Hilbert space Ti. we will label the generalized states by iply . In what 
follows, we consider the Dirac operator H = S in H = L2{^)'^. Finally, notice that 
Theorem 14.31 implies that the operator S defined on C^(f2)^ is essentially self-adjoint 
and has a unique self-adjoint extension H = S. 

In preparation of the proof for the completeness of the Chandrasekhar separation of 
variables ()2.32al -M for the Dirac equation in the Kerr-Newman geometry, we show first 
that the angular eigenfunctions form a complete orthonormal basis in L2(S'^)^. 

Lemma 5.1 For every keZ, j e Z\{0} the set {Y^^{^,ip)} with 

(5.2) 

is a complete orthonormal basis for L2{S'^Y . 
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Proof. Concerning the orthonormahty, we show that for every k' , k E and j' , 
j G Z\{0} the relation 

holds where (-I-) 52 denotes the scalar product on S"^. For simplicity in notation we omit 
to write the subscript uj attached to the angular eigenfunctions when there is no risk of 
confusion. With the help of ()5.2j) a direct computation gives 



52 



S^^, I sin 7? S^\'d)S''^ {'d). 



:<k j 



It suffices now to consider the case k = k . In this case since S''^ and S''^ are both 
eigenfunctions of the same differential operator U defined in (jS.lj) . they are orthogonal. 
Concerning the completeness, we can apply the projection theorem (e.g. Reed and 
Simon Theorem II. 3 vol I), i.e. we show that in L2{S^Y ^^e only element orthogonal to 
our orthonormal basis is the zero vector. Without loss of generality let us consider 




je<^+3)^ with G Co°"(0,7r) for 



1,2. 



After the variable transformation x = cost? we get 



52 



dx {ipi{x)S+{x) + ip2{x)S-{x)) 



Moreover, the following estimate holds 



52 



< 



dx (|<^i(a;)| \S+{x)\ + \<^2{x)\ \S-{x)\) 



1 



< llv^illt 

»i 



dx \S<(x) 



-1 



m\\L2 



dx \S^(x) 



<d dx (is+(j)r + \s-ix)\') = d 



where in the second line we used Holder inequality, in the third line the orthonormahty 
condition for S±{'d) and d := maxdlty^ilH^, ||¥'2||i2}- Clearly, | (^|F'^-^)52 1 < if and only 
if d = 0. Since d = implies v^i = v^2 = 0, C^{0,7i) is dense in L2(0,7r) and the scalar 
product (■|-)52 is continuous, the proof is completed. □ 

We state now the theorem on the completeness of the Chandrasekhar ansatz. In what 
follows let cr{H) C M denote the spectrum of the self-adjoint Hamiltonian operator H. 

Theorem 5.2 For every ip G C^{il)^ 



where the scalar product (-I-) is given by \4-b\l and x = {u,i!},ip). 



"^(■f^) jGZ\{0} keZ 
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Proof. We show first that it is possible to construct isometric operators 



W, 



. poo 



such that 



[u] 



u 



R^Uu) 



I R'JAu)Y^'A^,v) \ 

r'ja^)y!:a^,v) 



r:',my:aav) 

\RtU^)Y^'Ui^,v) J 



with some function A to be determined and given by ()5.2|) . Indeed, since the 
angular eigenfunctions Y^^ are normalized, we have 



2 

W \\L2{ 



du / d{cosi!}) / difip^ip^ 



2tt 



kj 



( p^^ 


ykj 

- -' UJ,~ 


\ 








Rt^ 






\Rt- 




) 



and by choosing A = (S) 4 it results that H-Rfi,-' || ^2(1^)2 = \\'^k,j{Rt^)\\'j^^(Q)4- By means 
of the isometric operators Wkj we can now introduce for every uj G cr{H) an auxiliary 
separable Hilbert space l){uj) as follows 



k,j& 



Moreover, the expansion theorem (e.g. Th.3.7 in Weidmann (1976)) implies that every 
element ■j/'^j in ^iyj) can be written as 

Notice that since the solutions b!^± of the radial system ()2.36p oscillate asymptotically 
for u — > +00 (see Section^), l)(a;) is not a subspace of L2(fi)'^ with respect to the spatial 
measure. Finally, let us recall that the direct integral of Hilbert spaces 



du f)(co') 



(5.5) 



a{H) 



is defined (see for example Yafaev Ch.l §5.1) as the Hilbert space of vector valued 
functions tp^u that take values in the auxiliary Hilbert spaces i){uj). By definition our 
Hilbert space Ti. will be decomposed into a direct integral ()5.5p if there is given a unitary 
mapping JF of 7i onto i^, but since the Hamiltonian H is self-adjoint, the spectral 
representation theorem (e.g. Theorem 7.18 in Weidmann (1976)) implies the existence 
of such an isomorphism and this completes the proof. □ 

Starting from the above result we get now more information about the spectral nature 
of the Hamiltonian H. In what follows we denote by aac{H) the absolutely continuous 
part of <j{H). 

Lemma 5.3 cr(iJ) is purely absolutely continuous and cr{H) = R. 
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Proof. We saw in Theorem 15 . 21 that every representative element ipi^ = {Til)){uj) of the 
element ip E Ti vn the decomposition ()5.5|) can be written in the form given by fl5.4|) . 
Hence, in order to investigate the spectrum of the Hamiltonian H we are allowed to focus 
our analysis on the radial system ()2.36p arising from the Dirac equation in the Kerr- 
Newman geometry when we make the ansatz ()2.H2al -M. As a consequence we just need 
to study the spectrum of the differential operator TZ associated to the formal differential 
system fl2.36|) after it is brought into the form of a Dirac system of ordinary differential 
equations (see Weidmann (1987)). Since this analysis has been already performed in 
Belgiorno and Martellini (1999) where it was shown that <y{TZ) = CadTZ) = K, we can 
conlcude from Theorem 15.21 that cr{H) = aac{H). 

Remark. A straightforward implication of the previous lemma is the absence of 
pure point spectrum for H i.e. app{H) = which means that no stationary states 
around a non-extreme Kerr- Newman black hole can exist. Regarding the non-existence 
of normalizable time-periodic solutions of the massive Dirac equation in the above 
mentioned framework see also Theorem 1.1 in Finster et al . (2000) and proposition 7.1 
in Hafner and Nicolas (2004) for the case of massless Dirac fields in the Kerr geometry. 
Hence, such black holes cannot form a new kind of atomic system with the charged 
black hole as its nucleus and around an electronic cloud. This possibility is instead left 
open in the case of an extreme Kerr- Newman black hole (see Schmid (2004)). 

In the next theorem we give the integral representation of the Dirac propagator in 
Kerr-Newman black hole manifolds. 

Theorem 5.4 For every ip G C^{Q)^ 

/+00 
iez\{o} k&z 

with scalar product (-I-) defined by l{4-6\) andipj'{x) as in / 1 5. ,9)) . 

Proof. Since H is self-adjoint. Theorem 7.37 in Weidmann (1976) guarantees that 
{U(t) = e^'^^lt G M} is a strongly continuous one-parameter unitary group with 
U{t)il) G C^{VlY for every t G M and every ip G C^(fi)^. The spectral theorem (e.g. 
Yafaev Ch.l §4-5) implies that 

/ + 00 
due'^\ct>^\i,^) (5.7) 
-oo 

where {(t)u)\ipLo) is given in terms of the resolvent of H by the following relation 

M.) = ^{ME{uj)i;) = ^ lim(0| [R{uj - te) - Rico + te^ (5.8) 

with E{uj) the spectral family associated to H. In addition Theorem 1.7 in Kato (Ch.lO 
§1) assures that {(puili^oj) is absolutely continuous in u, while the unicity of the spectral 
family E{uj) and the existence of a spectral family F{uj) = J^E{uj)J^^^ on follow 
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directly from Theorem 7.15 in Weidmann (1976). Notice that T is the unitary mapping 
already defined in Theorem 15.21 By means of equation (3) in Yafaev Ch.l §4.2 



30 . 

R{z) = I dE{u) 



we get 



1 1 /"+°° e 

[R^^ _ ,e) - R{uj + ze)] = - / ——^dE{u 



2-nV 

The above integrand is bounded and integrable. Hence, we may apply Lebesgue 
dominated convergence theorem to take the limit e — inside the integral sign. Since 

-lim- — - = 6{uj - u), 

it results that 

— lim[R{uJ - ie) - R{uj + ie)] = Id. (5.9) 
Finally, ^ follows from (Q, ^ and (Q. □ 

Remark. An integral representation for the Dirac propagator in the Kerr-Newman 
metric can also be found in Finster et al . (2003). In the derivation of the propagator 
they do not relate directly the spectral measure to the solutions of the radial and angular 
systems fl2.36p and ()2.37p but they consider first the Dirac equation in a finite volume 
[mi,'U2] X S"^ where certain Dirichlet boundary conditions on the spinors are introduced. 
With respect to our notation for the radial solutions these boundary conditions reduce 
to the following boundary conditions for the radial functions 

R4u2) = R4u2), R+{ui) = R-{ui). (5.10) 

The reason why they do that is that the spectrum of the Hamiltonian is purely discrete 
in a finite volume and hence, they are allowed to use the discrete version of the spectral 
theorem. Then they show that the spectral gaps for the Hamiltonian can be made 
arbitrary small in the infinite volume limit Ui ^ — oo and U2 —>■ +oo and this finally 
leads to the integral representation for the propagator. At this point notice that in such 
an infinite volume limit ()5.10|) is in contradiction with the asymptotics of the radial 
solutions in the next section. 



6. The asymptotic behaviour of the radial solutions 

In this section we give locally uniformly in u the asymptotics of the solutions of the 
radial system 



y/AV_ -imer{u)~x \ f R^{u) 
irrieriu) - A ] \ R+{u) 



0, (6.1) 



V± = — T K{u) = ij{r\u) + a^) - eQr{u) + ik + -\ a. 
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for u — > +00 and for u —00, respectively. Notice that the above system has been 
obtained from ()2.36|) by expressing the derivatives with respect to the radial variable r in 
terms of the tortoise coordinate u introduced in Section IHl Concerning the asymptotics 
for u +00, we treat the case |ti;| < me and |ti;| > mg, separately. Moreover, for ease 
in notation we will write r instead of r{u) when no risk of confusion arises. 

Lemma 6.1 Every non trivial solution R of 46'. 1\] for |c<j| > behaves asymptotically 
for u +00 like 

'coshe sinhe A / e-^*(")[/^ + 
sinhe coshe ) I e+^*(")[/^ + 



R{u) 




where f^ = c-te'''** for some uq > and the scalars scalars c_, c+ are such that 
|c_ p + |c+p 7^ 0. Moreover, 

1 fuj + me\ ^, . , Mml-ujeQ ,— ^ 

B = - log , '^[u) = Ku -\ log u, K = €[uj) Wuj^ — m'i 

with sign function e{ijj) such that e{ijj) = +1 if 00 > nif. and e{uj) = — 1 if uj < nif.. For 
\uj\ < rrie and uq > the asymptotic behaviour is 



R{u) 



with f^ = c'lpe^**-""^ for some scalars c_ and c+ such that |cLp + |c+p 7^ and 

'^[u) = ^/m^ — uo^u - log u. 

^Jm\ — up' 

Proof. Let us bring ()6.ip into the following form 

f = V(u)R(u), V(u) = ( ,^(;) I (6.2) 



with 



eQr {k + \)a \fK{X + irrier) 



+ + + d 

Since V{u) converges for u — > +00 to a constant matrix 
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—too inie 
—imp iuj 



we can decompose V{u) as follows 



V{u) = V^ + — + B{u) 
u 



with 



ieQ X - iMme \ \f^( \\ ^ 
X + iMme -leQ j' l^^^^^l " 
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for some positive constant C . Let us now define a matrix function %{u) as follows 



= V^o + — , = lim 



Notice that 



for M > Mo > the matrix functions and B{u) are continuously differentiable 
and 



+00 



du^Qi {u)\ < 00, 



+00 



(iu|S(r)| < cxD, 



• for leu I > me the eigenvalues Ai = —i^uj^ — m\ and A2 = +1^00"^ — ml of are 
simple. 

Let be an eigenvector of \4o belonging to the eigenvalue Aj and let \i{u) denote the 
eigenvalue of 2t(u) which converges to Aj as u ^ +00, then Theorem 11 Ch.IV in Coppel 
implies that the system ()(i.2|l has solutions 

u 



Ri{u] 

A simple calculation gives 



f" ds\i(s) 
gJ UQ ivy 



AiH 



6 + 

Mm2 - ujeQ 1 



1,2. 



A2(m) = -Ai(m). 



with K := H-a/cj^ — mg, while the eigenvectors of \4o are 



1 



,^1 = (cosh 9, sinh 0)"^ ^2 = (sinh 9, cosh 0)"^ , ® ~ ^ ^"-"S 



^uo — me 

Hence, for some scalars c_ and c+ every non trivial solution R{u), i.e. both c_ and c+ 
can not be zero at the same time, has asymptotically for u +00 the form 



R{u) = Gitu) 



cosh sinh B 
sinh cosh B 



where 



$(.):=..+ '^"^--"^^ log., /-:=c_e-(-), 



-j*{M0) 



The case \uj\ < m^. can be treated analogously. □ 



Remark. Unlike Lemma 3.5 in Finster et al . (2003) our method allows us to compute 
explicitly the expressions for the terms f '. 



00 
± • 



We give now an alternative proof of Lemma 3.1 in Finster et al .(2003), concerning the 
asymptotics of the radial solutions close to the event horizon. 

Lemma 6.2 Every non trivial solution R of \6.1\) behaves asymptotically for u ^ —00 
like 



R{u) 



[f'y + 0(6^^")] 

[fl + 0{e'^)] 
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with = ci°^e^^^°"i such that + \ff^\^ ^ and 

a f A; + — eOri ri — rn 



•where k+ zs i/ie surface gravity at the black hole horizon. 

Proof. Let us bring again ()6.1|) into the form ()6.2|) . Taking into account the definition 
of the tortoise coordinate u introduced in Section El it is easy to verify that 

u{r) = log {r -ri) + 0{r -ri), < k+ = J\ ^° (6.3) 

where k+ denotes the surface gravity at the black hole horizon. Taking into account 
that the potential V{u) entering in ()fi.2|) converges for u — oo to a constant matrix 

T/ r T/r ^ ( -'^0 ^ \ n J. g (fc + ^) - eQr, 
Vo = hm V{u) = , = ^ + — 5 ' 

fl6.3|) suggests an asymptotic expansion of V{u) in powers of 6^'^+". A short computation 
gives 

V{u) = Vo + O{{e^^+'')^). 

and applying Theorem 1 Ch.IV in Coppel it results that for u — oo the system ()6.2p 
has solutions 





for some ui G M such that u < ui < 0. Hence, for some scalars c?^ and c^"* every non 
trivial solution R{u) has asymptotically for u —oo the form 

R(u) - I ""^"""t^- ' + ^ f(°) - c(°)e^^^o- 

This completes the proof. □ 



Remark. Unlike Finster et al . (2003) by means of our method we get analytical 
expressions for the terms f^. Moreover we can show that their constant d entering in 
the asymptotics for the radial solutions at the event horizon can be expressed in terms 
of the surface gravity at the black hole horizon. This implies that d does not depend 
on ijj since depends only on the parameters a, Q and M. Thus the assertion in the 
statement of Lemma 3.1 in Finster et al . (2003) that d can be chosen locally uniformly 
in u is somewhat misleading. 
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7. The Dirac equation in oblate spheroidal coordinates 



In Section 12.11 we saw that by setting M = Q = the Kerr- Newman metric goes over 
to the Minkowski metric in oblate spheroidal coordinates ()2.7|) since by means of the 
coordinate transformation ()2.8|1 the expression ()2.7j) can be brought into the form of the 
Minkowski metric in cartesian coordinates. Concerning the computation of the Dirac 
equation in the Minkowski space-time in oblate spheroidal coordinates (OSC), we can 
again use the Carter tetrad ()2.9al -c) with A replaced now by r^ + a^ and proceed exactly 
as we did in Section 12.11 to end up with 



Oosc'^osc 





5_ 





-me 
5+ 



-me 








OSC 



0. 



(7.1) 



with 



a± = ±i 



a± = ±i 



(P± + /(r,^9)), 



+ 



(l)± + /(r,^?)), (3, 



V - — 



d ad 
dt + dip 



d 



(£± + ^(r,^)) 



d 



C± = — H — cot -(9 =F i a sin 'd— + esc '(^t— 



d 



M 



1 



+ a^ 



dt 
ia sin d 



dip 



2{r + ia cost?) 



r + ia cos'd ^ 

The Dirac equation 1)7.11) can be replaced by a modified but equivalent equation 

y^osci'osc = (^S''^OoscS)'$osc = 0. 
With the help of Corollarv 12.31 we can decompose Wosc as follows 



^{OSC) 



( 



{OSC) 



\ 



irneT 


Vr2 + a'^V. 






—irneT 






—irrier 








+ a2l?_ 


irneT 



J 



+ a^V+ 

where W(t,^,^) is given by ()2.19j) . Notice that the angular operator for the 

Dirac equation in the Kerr-Newman metric does not change when we go over to the 
Minkowski metric in OSC since it does not depend on the black hole parameters M and 
Q. If we now apply the Chandrasekhar ansatz as we did in Section 12.2) the equation 
yV^^^'^^'tposc = separetes into the angular system ()2.37|) and into the radial system 

j^OSC 
f>OSC 



y/7^Ta^V_ 



'inier 



inier — A y/r"^ + a^T) 




0. 



(7.2) 
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Concerning the construction of a symmetry operator for (|7.ip we can proceed as in 
Section flT^ Analogously to the procedure in Section HI we bring W^'^^'^^^osc = into 
the form of a Schrodinger equation idti^osc = Hosci^osc- Taking into account Sectional 
and that the Minkowski space-time in OSC can be continued analytically through the 
disc TZ defined in Corollary EHl it can be verified that the Hamiltonian Hqsc is formally 
self-adjoint with respect to the positive scalar product 



/ + 00 l*\ />27T 

dr / rf(cos^9) / dLp J " ^ i^osc4>osc- (7.3) 
-oo J-1 Jo y r + a 

In what follows for simplicity in notation we write ip instead of iposc- We consider now 
the free Dirac Hamiltonian Hqsc in the Hilbert space 

n = L2{nY := La ^ ^ dr d{cos^) d^ | , 1] = M x 5". 
From Theorem 1.1 in Thaller the operator Hqsc defined on 

4 



c^{nf = I n, \l ^ir^ d{cos^) dip 

is essentially self-adjoint and has a unique self-adjoint extension on the Sobolev space 
W^'\ny = W^'^ (n, \[^^-^ dr d{cos§) d^ 

In addition, the spectrum of Hqsc denoted by (Tosc, is purely absolutely continuous and 
given by (Tosc '■= (— oo,me] U [me,+oo). We state now the following results without 
proof since the proofs of Lemma 17.11 -3 are analogous, respectively, to the proofs of 
Theorem 15.21 Theorem 15.41 and Lemma f6. 11 

Lemma 7.1 For every ip e C(f 

J cTnf^n 2^-T?\ rr,i u^r^ 



where x = {r,i),(p), the scalar product (-I-) is given by (|7.5| ) and R^^± are the radial 



''^osc jez\{o} feez 



solutions of the system |yj 
Lemma 7.2 For every ip G C'^{VlY 

J ('OSC jez\{o} fcez 

with scalar product {■\-) defined by l{7.^ and tpj^lx) as in Lemma \7. 1\ 
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Lemma 7.3 Every non trivial solution R of (|7.ij| ) behaves asymptotically for r +00 
like 



R{r) 




coshe sinhe \ / e-*'^^[/~ + C»(r-i)] 
sinhe coshe ) I e+'''''[f^ + 0{r-^)] 



where f^ = c±e^*''''" for some tq > and scalars c_ and c+ such that |c_p + |c+p 7^ 0. 
Moreover, k, and are defined as in Lemma \6.1\ 
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